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The nonlinear Schrodinger equation (NLSE) in the presence of disorder is considered. The dy- 
namics of an initially localized wave packet is studied. A subdiffusive spreading of the wave packet 
is explained in the framework of a continuous time random walk. A probabilistic description of 
subdiffusion is suggested and a transport exponent of subdiffusion is obtained to be 2/5. 
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In this work the dynamics of an initially localized wave 
function is considered in the framework of the nonlinear 
Schrodinger equation in the presence of disorder. The 
wave function is governed by the following equation 



(1) 



where /3 is a nonlinearity parameter. The variables are 
chosen in dimensionless units and the Planck constant is 
h = \. The random potential V = V{x), x € {—oo, +00) 
is such that for the linear case {(3 — 0) the Anderson 
localization takes place P, Q , and the system is described 
by the exponentially localized Anderson modes (AM)s 
^uik = '^k{x), where '^^^^{x) are real functions and the 
eigenspectrum ujk is discrete and dense [s*]. The problem 
in question is an evolution of an initially localized wave 
function = 0) = "00 (2;) • 

This problem is relevant to experiments in nonlinear 
optics, for example disordered photonic lattices 0, [H, 
where Anderson localization was found in the presence 
of nonlinear effects, as well as to experiments on Bose- 
Einstein Condensates in disordered optical lattices [1, 0, 
H, y, [lOi [ll| . A discrete analog of Eg. (P) is extensively 
studied numerically [H, [H, [Tj, [l^, and a subdiffusive 
spreading of the initially localized wave packet has been 
observed, such that {x'^{t)) = J \ip{t)\'^x'^dx ~ t", where 
a transport exponent a was found in the range 0.3 0.4 

MM- 

A subdiffusive spreading of the wave packet was also 
obtained analytically 16] in the limit of the large times 
asymptotic dynamics of the tails of the wave packet, 
where the transport exponent was found to be a < 1. 
In that case the dynamics of the wave packet has been 
approximated by the fractional Fokker-Planck equation 
(FFPE) due to the arguments of a so-called continuous 
time random walk (CTRW). 

The concept of the CTRW was originally developed 
for mean first passage time in a random walk on a lattice 
with further application to a semiconductor electronic 
motion [ITj]. The mathematical apparatus of the frac- 
tional CTRW is well established for many applications 
in physics, see e.g., [H, [H, H [HI HI . 

The primary purpose of the present analysis is to de- 
velop the physical mechanism of the wave packet spread- 
ing and to obtain the transport exponent a. The anal- 
ysis is based on mapping the nonlinear Eq. ([T|) to the 



linear Liouville equation for the probability amplitude 
\ip{x,t)\'^, where the transition elements in the Liouville 
operator are determined by the nonlinear term in Eq. ([1]). 
Therefore we proceed to developing the analysis further 
by application of the CTRW approach to the correspond- 
ing Liouville equation (23j . 

First we obtain the linear Liouville equation for 
\'ip{x,t)\'^ [Hi- Projecting Eq. ([T]) on the basis of the 
AMs 

^P{x, t)=J2 Cu. (i)*-. (x) ^ Ckmuix) , (2) 



we obtain a system of equations for coefficients of the 
expansion Ck 

idtCu^uJkCk+P J2 Al'Xci.CkM,. (3) 



Here A(k) = AjJ \ is an overlapping integral of the four 
AMs: 



4;1 = / ■^k{x)^kAx)^M-^k.Ax)dx. (4) 



The initial conditions for the system of Eqs. ([3]) are such 
that ^^0(2;) ~ '^k^k'^k{x). Equations ^ correspond 
to a system of interacting nonlinear oscillators with the 
Hamiltonian 

H = Y^UkClCk + f3YAllfC*,^c:^Ck,Ck, . (5) 



Therefore, Eqs. ^ are produced by the Poisson brackets 
{H, . . .}pB by means of the Liouville operator 



L=-{H,...}pB 



f dH d dH d 



I — \dCl dCk dCk dC* 



(6) 

Since LH = and H{{C,C*}) = H{{a,a*}), we obtain 
that the Liouville operator is an operator function of the 
initial values: 



L 



dH d dH d 

da^. dak dak da^. 



(7) 



and corresponds to an infinite system of linear equa- 
tions dtC — LC, where C — C({afc,a^}) = 



2 



(. . . , Cfc_i, Cfc, Cfe+i, . . .) is an infinite vector. Thus, the 
Liouville operator reads 



d 



da, 



^4 



(8) 



where c.c. denotes the complex conjugation. The Liou- 
viUe equation is valid for any functions of the initial con- 
ditions {ak^ a^}. In particular, introducing the function 
Fk,k'{t) — Ck{t) ■ Cl,{t), one has the Liouville equation: 



dtFk,k' {t) = LFkM' (t) , Fk,k' {t = 0) = F, 



(0) 
k.k' 



akak' 



Therefore, the probability amplitude 

V{x, t) = = Fk,k' mk{x)-^k' {x) , 



k,k' 



as a function of the initial conditions, satisfies the Liou- 
ville equation as well: 



dtV = LV . 
Here the initial condition is 



(9) 



V{x, t^Q)^ro{x)=J2 Ft%'^k{x)^k' {x) . (10) 



k,k' 



Finally, we obtain that the NLSE ([T]) is replaced by the 
Liouville equation ([5]), which is the linear equation with 
a formal solution in the exponential form 



Vix, t) = e^*Vo{x) = J2 ^k{x)^k' (x) J2 [^^" 

k,k' n=0 



F, 



(0) 



k,k' ■ 
. 

In what follows we consider the dynamics of the initial 
wave packet Voix) in the framework of the probabilistic 
approach, where the dynamics of the wave packet is con- 
sidered as the CTRW. Since the dynamics of the prob- 
ability distribution function (PDF) V{x,t) is governed 
by the same Liouville operator in Eqs. ^ and (fTTjl . 
the overlapping integrals A(k) play the dominant role in 
the wave packet spreading. As follows from Eq. pT|) . the 
overlapping integrals determine the spread of the initially 
localized wave packet Vq (x) over all the Anderson modes 
as transitions from one mode to another. Since all states 
are localized, these transitions between states determine 
the transitions (or jumps) in the coordinate space as well. 

These transitions due to the overlapping integrals can 
be considered in the framework of the formal probabilistic 
integral equation Than the formal solution of Eq. 
(Ilip can have a form of the master equation 



rix,t) = 



g{x,t;x't')V{x',t')dx'dt' , (12) 



where the Green function g{x,t; x't') can be determined 
from the analysis of the overlapping integrals based on 
the probabihstic approach. [23 |. 

According to the values of the overlapping integrals, 
we divide the transitions between the localized states into 
two main groups. The first one corresponds to the ex- 
ponentially small overlapping integrals and the second 
one corresponds to the strong overlapping between four 
AMs when the overlapping integrals are of the order of 
1. In the case of strong overlapping, the AMs form clus- 
ters, where the wave functions have the same averaged 
coordinates for each cluster. Consequently, all transi- 
tions inside one cluster do not lead to any appreciable 
differences in the coordinate space, and we regard these 
transitions as trapping of the wave packet, or a particle, 
inside this cluster. Contrary to that, transitions due to 
the exponentially small overlapping integrals between the 
AMs belonging to different clusters lead to a change of 
the space coordinates that can be accounted for. We call 
these processes jumps. Contributions of trappings and 
jumps to the wave packet spreading described by Eqs. 
© and pT|) are different, and correspond to different 
probabilistic interpretations. 

In the sequel we follow the CTRW approach ^ and 
paraphrase it from [Tol . [2ll | in a form suitable for the 
present analysis. First, we consider a process of jumps. 
Let Pn{x) be the PDF of being at x after n jumps. It is 
reasonable to assume that the transitions between differ- 
ent states are independent of each other; therefore, the 
jumps are independent and obey the Markov property 



Pn+l{x) = / Pn{x')p{x - x')dx' , 



(13) 



where Po{x) = Vo{x) and p{x) is the PDF of a jump 
determined by the overlapping integrals as p{x) — 
^ exp(— ^|a;|)/2, and ^ — 1/R is an inverse localization 
length. 

The trapping is associated with clusters with effective 
lengths A. Due to the exponential localization, these 
values are distributed by the exponential law Pel (A) = 
A^^ exp (— A/Aq), where Aq is the mean length of a clus- 
ter. The effective lengths are determined by overlapping 
integrals; therefore, the minimum length of the cluster is 
A = _R, while the maximum one is A = AR. As a result, 
one obtains Aq — 5i?/2. Now the probability that a par- 
ticle exits this cluster and jumps to another one is of the 
order of ~ exp {—A/R). This value is also proportional 
to the inverse waiting time, t ~ exp (A/R). The proba- 
bility to find the waiting time in the interval (t, t + dt) is 
equal to the probability to find the corresponding trap- 
ping length in the interval (A, A -I- dA). Therefore the 
PDF of the waiting times is 



(14) 



where a = 2/5 and r is a time scale related to the trap- 
ping (25} . It follows that the mean waiting time is infi- 
nite. Taking into account that the waiting time PDF is 



3 



normalized we have 
wit) 



Wo 



l+a 



(15) 



such that w{t)dt = 1 and tw{t)dt — oo. Here 

2asin(7r/2Q) • t j_' ^ i. 

■"^o = ^(i+a)/2e. IS a normahzation constant. 

Now let us consider the PDF w{t) taking into account 
the dynamics of the jumps. Again, since transitions be- 
tween different states are statistically independent, wait- 
ing times for different jumps are statistically independent 
as well. Therefore, indexing the waiting time PDF by 
the jump number, we define that w„(t) is the probab ility 
density that nth jump occurs at time t (see e.g., [2l|, 
p. 42). Due to the reasonable assumption that jumps are 
independent transitions, we also introduce the Markov 
property for w„(i), which reads 



(t) = / Wn{t')w{t - t')dt' , 



(16) 



where wi{t) = w(t). Now we introduce the PDF 
P{x,t) — X)n arriving at coordinate x at 

time t. From Eqs. (fTS]) and (fTH]) we introduce the equa- 
tion [lii 

pix-x') / w{t-t')P{X' ,t')dx'dt' 

-oo Jo 

+ Vo{x)6{t), (17) 

which relates the PDF P(x, t) of just having arrived at 
position X at time t to the PDF P{x',t') of just arriving 
at x' at time t'. The last term in Eq. p7|) is the initial 
condition. Thus the PDF V{x,t) of being at position x at 
time t is given by arrival at x' at time t' and not having 
moved after this event, namely 



V{x,t) 



P{x,t')W{t-t')dt' , 



(18) 



where W{t) = 1 — J^w{t')dt' denotes the probability of 
no jump during the time interval (0,t). Performing the 
Fourier transform p{k) — J-p{x) and the Laplace trans- 
form wis) = Cw{t), we obtain the MontroU- Weiss equa- 
tion [17] 



Vik, s) = TCP 



l-w{s) Poik) 



1 — p(k)w{s) 



(19) 



This expression determines the master equation (|12p 
and establishes a relation between the Green function 
and the overlapping integrals in the CTRW form: Q — 
p{x — x')w{t — t'). Eq. p9|) can be simplified for the long 
time s ^ 1 and the large scale fc ^ 1 asymtotics that 
corresponds to the diffusion limit (fc, s) — > (0,0). Tak- 
ing into account the Fourier p{k) and the Laplace w{s) 
images in Eq. p^ : 



w{s) 



1 



1 



1 + (sr)" 



1 - R^k^ , 
. 1 - {stT 



(20) 



we obtain for the PDF in the Fourier-Laplace domain 
(see also [ll|) 



P{k,s) = 



Po{k)/i 



l+DaS~°'k^ 



(21) 



where Da = P? /t^ is a generalized diffusion coefficient. 
Using the Laplace transform of the fractional integration 



J{r)di 



(t-r)i- 

one obtains the fractional integral equation 

P{x,t) - Poix) = di''DadlP{x,t) . (22) 

Differentiating this equation with respect to time, one 
obtains that the CTRW is described by the fractional 
Fokker-Planck equation (FFPE) ^ 



dtP{x, t) - Dadr"diP{x, t)^0 



(23) 



where is a designation of the Riemann-Liouvillc frac- 
tional derivative 



dt 



1 



f{r)dT 
r(l - ly) dt Jo (t - Tf 



where Q < v < 1 . From Eq. one obtains for the sec- 
ond moment (x^(i)) = J_^x^P{x,t)dx the differential 
equation: 



2Dat° 



Via) 



Here T{z) is the gamma function, x{t = 0) = 0, and 
we use the following property of the fractional deriva- 
tive di'l = t-^/Til - v). Therefore, Eq. ^ describes 
subdiffusion 



{x\t)) = 



2Dat°' 

T{l + a) 



(24) 



with the transport exponent a — 0.4. In the recent nu- 
merical studies of the discrete NLSE [11, HI] the exponent 
a was found in the range 0.3 0.4. 

This consideration can be extended on the wave packet 
spreading in the framework of the generalized nonlinear 
Schrodinger equation 



(25) 



where n > 1 and the Hamiltonian TYq has the pure point 
spectrum with the localized eigenf unctions: Ho^fc = 
Sk'^k- For example, the Hamiltonian describes Wannicr- 
Stark localization [23], and the discrete counterpart of 
Eq. ((25)) with n = 1 corresponds to delocalization in a 
nonlinear Stark ladder [2^ [2^ . Repeating probabilistic 
consideration of the CTRW based on the overlapping in- 
tegrals A(k) of 2ri -t- 2 eigenfunctions ^'^(a;), one obtains 



4 



that Eq. I|25p describes subdifFusion with the transport 
exponent 



For n ~ 1 this expression coincides with the numerical 
result of Ref. [H- 

The nonlinear Schrodinger equation in the presence of 
disorder is considered, and the dynamics of an initially 
localized wave packet is studied. A subdiffusive spread- 
ing of the wave packed is explained in the framework of a 
continuous time random walk. It is shown that the sub- 
diffusive spreading of the initially localized wave packet 
is due to the transitions between those AMs which are 
strongly overlapped. This overlapping leads to the clus- 
tering with an effective potential well and, correspond- 
ingly, to effective trapping of the wave packet inside each 
cluster by the potential. Therefore, the dynamics of the 
wave packet corresponds to the CTRW, where the basic 
mechanism of subdiffusion is an entrapping of the wave 
packet with delay, or waiting, times distributed by the 
power law w{t) ^ < a < 1. The trapping 



mechanism also determines the transport exponent a, 
which is due to the number of AMs contributed to the 
overlapping integrals according to Eq. ([26|) . One should 
recognize that the PDF of the CTRW in Eq. ^ is as- 
sociated with the true quantum distribution of the NLSE 
by virtue of Eq. ^ with the formal solution of Eq. 
These expressions are an important intermediate result 
which establishes a relation between the NLSE and the 
probabilistic description. Note, that the PDF Vix, t) in 
Eq. ^ is the exact distribution, and it corresponds to 
the CTRW in the long time and the large scale asymp- 
totics of Eq. (PO)) . 

It is worth noting, in conclusion, that the linear prop- 
erty of the Liouville equation ^ and its formal solution 
of Eq. (jlip are important for the probabilistic approach. 
This linear property was used for the CTRW approach, 
and the MontroU- Weiss Eq. determines the Green 

function in Eq. 
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